m 
m 



A stochastic model for tumor growth with immunization 

Thomas Bose and StefFen Trimper 
Institute of Physics, Martin-Luther-University, D-06099 Halle Germant 

(Dated: September 1, 2008) 

■ Abstract 

o . 

(•~^ . We analyze a stochastic model for tumor cell growth with both multiplicative and additive colored 

<N , 

^jr^l noise as well as a non-zero cross-correlations in between. Whereas the death rate within the logistic 

< 

• be stochastically changed due to biological motivated growth processes leading to a multiplicative 

internal noise. Moreover, the system is subjected to an external additive noise which mimics the 
influence of the environment of the tumor. The stationary probability distribution Ps is derived 
Q ■ depending on the finite correlation time, the immunization rate and the strength of the cross- 

,jp . correlation. Pg offers a maximum which becomes more pronounced for increasing immunization 

rate. The mean-first-passage time is also calculated in order to find out under which conditions the 
tumor can suffer extinction. Its characteristics is again controlled by the degree of immunization 
. and the strength of the cross-correlation. The behavior observed can be interpreted in terms of 



model is altered by a deterministic term characterizing immunization, the birth rate is assumed to 



the three state model of a tumor population. 
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I. INTRODUCTION 



A fundamental aspect of all biological systems is the understanding of emergence of coop- 
erative behavior. The competitive interaction among different growth and death processes 
and the inclusion of external mechanism are widely believed to influence the global properties 

n 

of such systems [l]. In this context much effort has been devoted to model the dynamics of 
competing population through a nonlinear set of rate equations such as proposed by Lotka 



and Volterra [2, y| or a broad variety of their variants as a stochastic model for ecosys- 

n n n 

tems j^, coexistence versus extinction [5| or special clustering in Lotka- Volterra model [6[. 
Prey-predator systems are likewise related to that kind of models, where recently also fluc- 
tuations and correlations are discussed [7|, l9j as well as instabilities with respect to spatial 
distributions [8]. The heuristic approach is based upon deterministic evolution equations. 
Otherwise, a population of proliferating cells is a stochastic dynamical system far from equi- 
librium Proteins and other molecules are produced and degraded permanently. Cells 
grow, divide and inherit their properties simultaneously to the next generation. To gain 
some more insight into the generic behavior of phenomena such as tumor cell growth, it is 
desirable to take into account both internal and external stochastic noises as well as spatial 
correlations. 

In the present paper we are interested in tumor growth which had been attracted attention 
over several decades. Mathematical modeling of the growth of a certain population is based 
on different equations where the logistic growth and the Gompertz law are the most popular 



deterministic models 



A more refined model was presented in [l2], however we argue that 



the solutions for the stationary probability distribution and the mean-first-passage time are 
not calculated correctly. The details and the corrections are given in our paper in Sec. Ill and 
IV. Nevertheless, the model in jl^] includes already both additive and multiplicative noise 
terms considered likewise in 13| . However, the stationary distribution function presented in 



that paper is also not correct as pointed out in and replied in [l5|], see also our results 
discussed below. The role of pure multiphcative noise may induce stochastic resonance, 
which appears in an anti-tumor system Q- In that work the deterministic forces are 
modified as it will be also discussed in the present paper. The mean first passage time of a 
tumor cell growth is altered by cross correlations of the noise, see already jl^]. Essentially 
for tumor modeling is also the inclusion of therapy elements as proposed in [l^. In our 
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model we analyze a special immunization term which enhances the the death rate. Another 
model is devoted to the spatiotemporal triggering infiltrating tumor growth. 

Our approach can be grouped into the permanent interest in a statistical modeling of 
growth model, where evolution equations of Langevin or Fokker-Planck-type play an decisive 



role [20|. In particular, the focus is concentrated on correlated colored noises [2l| in the 



form of multiplicative noise 2^ and additive noise 23|] . A similar approach is also applied 



m 



24| for the BernouUi-Malthus-Verhulst model. In the context of population dynamics 
different aspects has been studied such as time delay effect s |25l| . a general classification 
scheme for_phenomenological universality in growth problems [26[ , extinction in birth-death- 
systems [27i], the complex population dynamics as competition between multiple-time-scale 



phenomena 28| and the the dissipative branching in population dynamics 29 1. 

The goal of our paper is inspired with the view to alter the models in such a manner 
that both immunization and correlated noise are included. Especially, we want to demon- 
strate that a finite correlation time and a nonzero immunization rate have an significant 
impact on the different steady states realized within the model. Additionally we analyze 
the interplay between an internal noise leading to a stochastic birth rate and an external 
noise. Furthermore, the mean-first passage time is calculated which enables us to analyze 
under which conditions, depending on the correlation time and the immunization rate, the 
tumor population can suffer extinction. Our paper is organized as follows: In Sec. II we 
define the Langevin equation with different multiplicative noises and their cross- correlation 
functions, the meaning of that is considered in detail. Then we introduce an immuniza- 
tion term the influence of which will be analyzed in the paper. Such an additional term 
leads to a significantly modified death rate. Based on the related Fokker-Planck equation 
the stationary probability distribution (SPD) is studied in Sec. III. The expression for the 
mean-first-passage time (MFPT) and its meaning is explained in Sec. IV. Further we discuss 
the relation of our results to real tumor growth. In Sec. V we finish with some conclusions. 

II. THE TUMOR MODEL 

In order to develop a statistical tumor growth model, we consider the general type of 
Langevin equation, that reads 

dx 

= fix) + giix)€iit) + (/2(x)e2(t), (1) 



where x{t) denotes the number of tumor cells at time t, f{x),gi{x) and g2{x) are deter- 
ministic functions and ei(if:) and e2(t) are colored noises with zero mean and colored cross- 
correlation. These statistical properties are given by (ei(t)) =0; (e2(t)) = and the 
corresponding correlation functions 



C{t - t') 



{e,{t)e^{t')) {e,{t)e,{t')) j \^ exp ^ exp 




T2 

Here, the elements of the correlation matrix Cij{t—t') are assumed to be symmetric Cij 
The quantities M and a are the noise intensities and ri and T2 are the correlation times 
of the autocorrelation functions Cu and C22- The parameters A and T3 characterize the 
strength of the cross-correlation function between ei(t) and £2{t) and the cross-correlation 
time, respectively. In our model we consider a modified logistic growth model with 



f(x) = ax — b(x)x , b(x) = &o + 4>ix) = ^0 H o (3) 

1 + x^ 

Here, the parameter a is the deterministic growth rate and 60 denotes the decay rate propor- 



tional to the inverses carrying capacity, respectively Ms deatl> rate is^altered by inclusion 

of a tumor- immunization interaction represented by the function <f>{x) [1], where the param- 
eter (3 designates the strength of the immunization. Under immunization the effective death 
rate b{x) is enhanced where the decay of the rate depends on the immunization strength 
(3. The behavior of the effective death rate is depicted in Fig. [H The tumor cell evolution 
is further coupled to internal and external noises denoted by ei{t) and e2{t), respectively. 
Whereas the death rate is systematically enhanced by immunization, modeled by the deter- 
ministic function (p{x), the effective birth rate should be influenced by the stochastic force 
ei(t). This leads to the assumption 

giix) = -x. (4) 

Furthermore, the effect of additive noise represented by e2(t) is incorporated into the system 

by 

g2{x) = l. (5) 

Notice that all parameters are dimensionless, so that the prefactors in the last equations 
could be set as unity. With regard to the discussions in Sec. IV let us introduce an effective 
potential V{x) according to the deterministic force f{x), that reads 

V{x) = - [ f{x)d{x). (6) 



b 1,4 




FIG. 1: Plot of the death rate 6 as a function of the cell number x when 6o = 1 is fixed: (3 takes 
(solid line), 0.2 (dotted line), 0.5 (dashdotted line) and 1 (dashed line). 



Evaluating ([6]) yields the following expression for V{x) : 

V{x) = -bo — -a x'^ + P {x — arctan x) . 



(7) 



The potential V[x) is presented in Fig. [2J The stationary points can be determined by 
setting /(x) = ~^^(3;) = 0. From here we discriminate between four extrema, from which 
only two of them are real in the parameter range considered. The remaining stationary 
points take complex values and will not be discussed furthermore. Thus, we have derived a 
potential with a minimum at xi = > and a maximum at X2 = 0. 



III. FOKKER-PLANCK-EQUATION 



A. Derivation of the stationary probability distribution (SPD) 



As a next step the Langevin Equation ([T]) is transformed into an equivalent Fokker-Planck 
equation (FPE) ^flflQ. To that aim let us consider x(t) as a random variable whose 
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FIG. 2: Plot of the effective potential y as a function of the cell number x where a = 1, bo = 1 
and [3 = 1. 

probability density function p{w,t) is a delta-function 

p{w,t) = 6{x{t) — w) . 

From here one can find the stochastic Liouville equation [s^ for the probability distribution 
function 

P{w,t) = {p{w,t)). (8) 

Here, P{w,t) is the density of the probability distribution function that the process x{t) 

takes the value w at time t. From this relation combined with Eqs. ([I]) - ([5]), one obtains 

the FPE in the form 

dP(w f) 8 (9^ 

= + ^,[B{w)Piw,t)]. (9) 

The explicit expressions for A{w) and B{w) are 

M \ /, , . 1 \ 2 Av^Ma 



A{w) = { a+ ] w - { bo + P- ^ w 

^ ' * l + arW V l + toV l + arg 



B{w) = w^-—-^ w + - . (10 

1 + ari l + ars 1 + ars ^ ^ 
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Notice there is a relation between the functions A{w), B{w) and f{w) of the form 

Aiw) = fiw) + ^^Biw). (11) 

The stationary probabihty distribution (SPD) of the system can be obtained from Eqs. 
( ITTl) ) and can be written as 



PJw) 



N 



'Biw] 



exp 



Biw') 



dw' 



where is the normahzation constant that is determined by 



PJw) dw = 1. 



(12) 



(13) 



Depending on the cross-correlation strength A one has to distinguish between different cases. 



The solution of the SPD for < A < 



PJw) 



l+OTS 



i/(l+ari){l+ar2) 



IS 



'B w) 



exp 



U{w) 
M 



(14) 



where we have introduced a generalized potential according to 



U{w) = h{w)-E\n [B{w)] 



Fi — F2 



VMa - A2 

Here, the following abbreviations are utilized: 

M 

M = , a = 

1 + ari 1 + aT2 

The non-universal exponent E reads 



arctan 



Mw - A 
VMa - A2 



a ~ AV Ma 

1 A = 

1 + aTs 



~ _ a(l + ari) Ha (1 + an)^ PMXVMa 

~ 2 ° yM 1 + arg ~ K{1 + org) 



with 



Further we use 



K 



M 



a 



1 + an 1 + aT2 



+ 



(15) 



(16) 



(17) 



(18) 



h{w) 



1 + ara 
(3Ma ( a 
K 



2A 



2 1 + ari 



M 



l + ara (1 + ars)^ 

2X'^Ma 



+ 



l + ar2)2 (l + aTi)(l + ar2) (1 + ars)^ 
M 



a 



1 + ari 1 + ar2 



AyMa r 21 1 

arctan w m 1 + to \ + y{ 

1 + aTs -' ' 



w) 



y{w) = 60(1 + an) w 



(19) 
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Let us remark, that by settine tt = To = and 3 = one obtains the corrected solution for 

n 

the SPD in [IJ] (equations (19)-(22)). In the present paper we assume that all correlation 
times take the same values, that is ti = r2 = T3 = r resulting in new expressions for the 
generalized potential denoted now as U{w). In case of the condition < A < 1 we get 



Ps{w) 



N 



where B(w) according to Eq. (fTOl) changes to 



exp 



U{w) 
M 



B{w) 



The potential takes the form 



M 



1 + ar 



2 ^ /\ \ ivi (y. 

W : W + 



a 



1 + ar 



1 + ar 



U{w) = h{w) - E\n [B{w)] 



F(l + ar) 



arctan 



Mw - \VMa 



with 



E 

Q 

F 

h{w) 
y{w) 



a(l + ar) 



a 




M^ + a'^ + 2Ma{2X^ - 1^ 



lAVMa - h^a {2X^ - l) 



AW^(l + aT) 



f3Ma {a + M(2A2 - 1)) 



Q 



(3M{1 + ar) 

Q 

60(1 + clt) w. 



({M - a) arctan w - \\fMa In [l + w'^]^ + y{w) 



(20) 



(21) 



(22) 



(23) 



We want to point out that setting [3 = and r = gives the right results for the SPD that 



is not correct in [l^. These solutions are in agreement with those mentioned in 1^. The 
second case of A = 1 has to be considered separately. The corresponding solution is 



Ps{w) 



N 



exp 



and the generalized potential reads 

U*{w) = h*{w) - E* In [B{w)] 
The non-universal exponent is written in the form 



U*{w) 
M 



F*(l + ar) 



y/Ma- M 



w 



E* 



a(l + ar) 



a 



(24) 



(25) 



(26) 



with 



Q* = {M + af 

17* h l3Ma{a + M) 

F = ay Ma — br,a 

Q* 



h*{w) = \t ((^ ~ arctanw - v^Afo In [l + w^] j + y{w) . (27) 



PM{1 + ar) 

The function y{w) remains unchanged and is given by Eq. fl23l) . In the following sections 

we only consider the first case and analyze the results for < A < 1, i.e. our further 
computations refer to the solutions given by Eqs. fl2U]) - (1251) . 



B. Properties of the SPD 

In this section we discuss the behavior of the stationary probability distribution (SPD) 
calculated analytically in the previous subsection. In Fig. [3] the SPD is represented as 
function of the tumor cell population w under different immunization rates (3. The SPD 
reveals a maximum indicating the most probable cell population. The maximum becomes 
the more pronounced the higher the immunization rate (3 is. The maximum is shifted to 
smaller tumor population with increasing rate (]. The SPD is influenced significantly by the 
cross-correlation characterized by the parameter A. The maximum is strongly enhanced by 
an increasing cross- correlation strength as shown in Fig. HI The SPD is also influenced by 
the correlation time r of the noises. The result is shown in Fig. [51 There appears already a 
maximum which is more articulated when the correlation time is enhanced. 



C. Biological interpretation 

The importance of an efficient immunization against tumor evolution is illustrated in 
Fig. [3l This efficacy depends on the competence of the immune system to detect the malig- 



nant cancer cells, and thus to initiate a power full immune response. In 3l|, l32] the 'Three 
E's of cancer immunoediting' are described, i.e. the tumor-immune interaction can result 
into three different phases: elimination, equilibrium and escape, whereas sooner or later 
the equilibrium phase offers a cross-over to the other phases. Adopting this concept to the 
behavior of the SPD the tumor elimination phase is the more probable and the escape phase 
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FIG. 3: Plot of the SPD Pg as a function of the ceh population w for fixed a = 0.5, 60 = 1.0, 
a = 0.3, M = 0.7, r = 0.5 and A = 0.5. The immunization /? varies from 0.0 (solid line), 0.5 
(dotted line), 1.0 (dashed line). 

is the more improbable the higher the immune coefficient (3 is, for further remarks compare 
also subsection IV. C. 

Now, we want to relate the internal noise ei and the external noise 62, introduced in 
Eq. ([T]), to real processes that occur in the tumor and its environment including the hallmarks 



of cancer 33|], and moreover to point out to the connections among each other. External 
noise is thought to be originated in the extracellular matrix embedding the tumor or it is a 
consequences of drug delivery from outside of the host. Additionally, external noise can be 
caused by thermal fluctuations. In contrast the internal noise is generated directly within 
the tumor system as a kind of self-organization, for instance by gene mutations resulting 
in a multitude of genetically different tumor cells within the same system. The process 
is based upon internal mechanisms inside the tumor without contact to its environment. 
Although the origins of both stochastic processes are different one should argue that there 
exists an interrelation among both ones. A measure for such a correlation is the strength 
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FIG. 4: Plot of the SPD Pg as a function of the ceh population w for fixed a = 0.5, 60 = 1.0, 
a = 0.3, M = 0.7, r = 0.5 and /3 = 0.5. The strength of the cross-correlation A takes 0.1 (solid 
line), 0.5 (dotted line), 0.9 (dashed line). 

of the cross-correlation denoted by A in Eq. ([2]) as well as the correlation time r. The 
expected coupling between external and internal noises can be understood as follows. The 
normal tissue adjacent to the malignant one produces anti-growth signals in order to avoid an 
uncontrolled growth. The tumor may respond with insensitivity with respect to these signals 
by alteration or down-regulation of the corresponding receptors. Furthermore, some tumor 
cells are able to develop self-sufficiency in generating growth signals. Another correlation 
concerns the nutrient supply. With a growing tumor tissue the competition is intensified 
regarding the nutrients between normal tissue and the nascent transformed cells. The tumor 
can sustain and induce angiogenesis via an 'angiogenic switch' from vascular quiescence in 
order to progress to a larger size. Another characteristic of tumor growth is the acquisition 
of a diversity of strategies to evade apoptotic signals that are emitted on the one hand by 
the tumor environment and on the other hand generated within the tumor cells. 

Therefore, the behavior of the SPD depending on the strength of the cross-correlation is 
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FIG. 5: Plot of the SPD Ps as a function of the cell number w when a = 0.5, 60 = 1-0, A = 0.5, 
M = 0.7, a = 0.3 and (3 = 0.5 are fixed: r takes 0.0 (solid line), 0.5 (dotted line), 1.0 (dashed 
line) . 

clearly shown in Fig. HI An increasing A is equated with an increasing ability of the tumor 
to compensate the external interferences via internal reactions described above. Thus, in 
case of strong correlations the tumor has an improved ability to reach the escape phase. 

In order to explain the dependence of our results on the correlation time let us remind 
that r is the correlation time of the cross- correlation as well as the correlation time of 
the auto-correlation functions of the additive (external) and multiplicative (internal) noise, 
respectively. Here we have assumed that the correlation time for both kind of noises is 
relevant on the same time scale r. Taking this into account the appearance of a finite 
correlation time leads to a higher probability of a certain tumor size but does not change 
the likeliest tumor size as presented in Fig. [51 

Notice that we attribute a random nature to the mechanisms of the tumor evolution 
because the details of the growth and decay processes differs from patient to patient. There- 
fore, tumor growth and the interplay with the environment can be regarded as a stochastic 
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process and is interpreted by introducing external and internal noises. 



IV. MEAN-FIRST-PASSAGE TIME (MFPT) 



A. Derivation of the MFPT 



In cancer treatment it is of interest whether a tumor that reached a certain size can suffers 
extinction by external or internal interferences, i.e. is it possible that the influences of the 
noises and the immune system introduced in the previous sections can cause extinction of the 
tumor. A further concern is the transition time between these two states: the lethal tumor 
size and the tumor free state, respectively. In order to describe these transient properties of 
;he system we apply the mean-first-passage time that is given by the following expression 



34 



35| 



W2 



T 



W1W2 



dw 



B(w)PJw) 



Ps{v)dv, 



(28) 



i.e. the transition from an initial point wi to an end point W2 is considered. We choose 
the stationary points of the effective potential (171), more specifically wi = Xs and W2 = 0, 
i.e. the MFPT of the system reaching the tumor free state is studied. Now we make use of 
an approximation scheme that is valid for small M and a in comparison with the potential 
carrier high [U{w2) — U{wi)] |20|, l36| what has already been applied to tumor models, e.g. 
~v\ . We derive an analytical expression for fl28l) . namely 



2n 



- W1W2 



exp 



- [^(0) - U{w,)] 



(29) 



v"io)v"iwi: 

where the double-prime denote the second derivation with respect to w. Inserting Eq. ([7]) 
and Eqs. fl22|) - fl23|) into the Eq. fl29l) leads to the final expression 



T - ^ 
i ^^exp<j^^ 

F(l + ar) 



+ 



arctan 



a 



arctan 



Xy/Ma - Mwi 



(30) 
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where 



R 




2Xy/M'a 



1 + ar 1 + ar 1 + ar 1 + ar 

= ^^^y" "'^^ (^(M - a) arctan wi - Av^Ma In [l + wi^] ^ + y{wi) 
y{wi) = boWi{l + ar). (31) 

Both constants, E and F, are still the same as those in Eq. (123!) . Notice, that applying our 
solutions obtained by Eqs. f|T5l) - f|T9l) into Eq. fl29|) . therefore substituting f/(0) and f/(wi) 
by f/(0) and f7(tyi), respectively, and setting ri = r2 = and (3 = yields the correction of 



the expression in Eq. (27) in 121]. 



B. Properties of the MFPT 

In this subsection we discuss the properties of our system and the behavior of the MFPT. 
In Fig. [6] we present the MFPT as function of the parameter M introduced in Eq. ([2]). This 
parameter M is a measure for both the auto- and the cross-correlation function between 
internal and external noise. As a feature there occurs a maximum indicating a long living 
cell population. The maximum is the more pronounced the lower the immunization rate is 
and it is shifted towards higher values of M. Increasing the rate (3 the MFTP is smaller and 
an extinction of the tumor population is more probable. In Fig. [7] the MFPT is represented 
depending on the parameter a according to Eq. ([2]). Here a characterizes the strength of 
the auto-correlation of the additive noise as well as the strength of cross-correlation. The 
increase of a leads to a decrease of the MFPT. This decay is very strong in case of a high 
immunization rate as expected. The direct influence of the immunization strength (] on the 
MFPT is shown in Fig. |8l There appears already a maximum which is shifted to higher values 
of P when the correlation time r is reduced. A similar behavior of the MFPT as function 
of P is also observed in dependence on the parameter A. A very instructive behavior can be 
observed in Fig. [9] where the MFPT is depicted as function of the immunization coupling 
(3 with variation of the global noise strength M. The maximum becomes more pronounced 
if the noise strength increases. A nearly linear behavior of the MFPT as function of the 
correlation time r is observed in Fig. [TOl The increase of the MFTP is weaker for a stronger 
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FIG. 6: Plot of the MFPT as a function of M when for fixed values o = 1.0, 6o = 1-0, r = 0.5, 
a = 0.1 and A = 0.8. The immunization strength /3 varies from 0.3 (solid line), 0.5 (dotted line), 
0.7 (dashed line). 

immunization rate j3. 

C. Biological aspects 

In this subsection the behavior of the MFPT is discussed with regard to biological aspects. 
Let us stress that a decrease of the MFPT is tantamount to an increase of the probability of 
the transition to the tumor free state. At first, we consider the influence of the multiplicative 
noise on the MFPT and its relation to the immune system, where the noise is originated 
from all the distinct processes described in subsection III.C. Here, we assume that the 
multiplicative noise is mainly determined by the gene mutations. Fig. [6] indicates that there 
exists an appropriate M leading to a high MFPT. Improving the effectiveness of the immune 
system leads to a reduction of the MFPT. Moreover, the tumor system requires more gene 
mutations in order to maximize the MFPT. More genetic alterations induce a deterioration 
of the ability of the immune system to identify tumor cells. But this mechanism is limited 
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FIG. 7: Plot of the MFPT as a function of a when a = 0.5, 60 = 1-0, r = 0.5, M = 0.8 and A = 0.5 
are fixed. The parameter (3 takes 0.1 (solid line), 0.5 (dotted line), 0.9 (dashed line). 

as it is visible by the descent of the curves in Fig. [HI As soon as the optimal value of the 
strength of the multiplicative noise is exceeded the MFPT decreases and consequently the 
ability of the self-organization is reduced. 

The influence of the external (additive) noise offers the expected behavior. All interfer- 
ences impair the living conditions of the tumor. Therefore, a growing parameter a leads to 
a decline of the MFPT and enhances the probability of the extinction of the cancer. 

Following the explanation made for the interpretation of Figs.[8]and [9]let us introduce the 



principle of immunoediting presented in 3l|, |32] . On the one hand the immune system is able 



to cause extinction of the tumor, otherwise it can facilitate tumor progression by sculpting 
the immunogenic phenotype of tumors. During this process immune-resistent variants of the 
tumor cells are able to survive and even more proliferate in order to develop a tumor tissue 
that can sustain further immune attacks. This behavior is displayed in Figs. [8] and [91 The 
increasing part of the curve is thought to be connected with the process of tumor sculpting 
which may end up in the tumor escape phase. The decreasing part of the MFPT is identified 
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P 

FIG. 8: Plot of the MFPT as a function of (3 for fixed values for a = 0.5, ho = 1.0, a = 0.6, M = 0.8 
and A = 0.5. The correlation time r varies from 0.1 (solid line), 0.5 (dotted line), 0.9 (dashed line). 

with immunosurveillance that leads to tumor ehmination. This process is the more hkely 
the bigger the immunization strength is. The effect of the strength of the muhiphcative 
noise M and the strength of the cross-correlation r are similar. Due to the fact that the 
increase of both parameters M and r and the biological mechanisms beyond promotes the 
tumor growth, one should expect a retardation of the transition to the elimination phase as 
shown in Figs. [8] and [H 

The correlation-time r also effects the MFPT. An increase of r leads to a slowing down 
of the transition between the different states of the tumor. The longer the correlation time 
r is the more probable are long living tumor populations. Consequently, a rising value of r 
simplifies the opportunity of the tumor to evade the immune system. 

V. CONCLUSIONS 

In this work we have proposed and analyzed a more refined model describing tumor cell 
growth. Starting from a logistic model we have modified the model in several directions. 
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FIG. 9: Plot of the MFPT as a function of (3 when a = 0.5, bo = 1.0, a = 0.6, A = 0.5 and r = 0.5 
are fixed. The noise strength M is 0.1 (sohd hne), 0.5 (dotted hne), 0.9 (dashed hne). 

The decay term is supplemented by a deterministic non-linear immunization term which 
enhances the death rate of the tumor. Furthermore, the birth rate as assumed to be stochas- 
tically distributed leading to a multiphcative noise. The occurrence of such a noise term is 
motivated by the underlying biological situation. Additionally, the system is subjected to 
an additive, external noise which is originated by the external conditions as the environment 
of the tumor. Both kinds of colored noises are correlated, i.e. there are autocorrelation and 
a cross-correlation functions with different strength. The resulting equation has the form 
of a Langevin equation which can be transformed into a Fokker-Planck equation. Using 
standard methods we find the steady state solutions which are discussed depending on the 
strength of the cross-correlation, the finite correlation time and the degree of immunization. 
The behavior of the stationary probability distribution is analyzed taking into account bi- 
ological aspects above all the three different states of the tumor: elimination, equilibration 
and escape phase. In particular, the SPD offers a maximum indicating the appearance of 
very probable states. This maximum becomes for instance the more pronounced the higher 
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FIG. 10: Plot of the MFPT as a function of r when a = 0.5, bo = 1.0, a = 0.3, M = 0.8 and 
A = 0.5 are fixed: /? takes 0.1 (sohd Hne), 0.5 (dotted hne), 0.9 (dashed hne). 

the immunization rate is. As a further quantity of interest we have studied the mean-first 
passage time which indicates when the tumor suffers extinction. The MFPT is hkewise 
calculated analytically and analyzed under consideration of biological aspects. The MFPT 
is influenced in a significant manner by the immunization and the cross-correlation as well 
as the finite correlation time of the underlying colored noises. The observed behavior can 
be understood in terms of the above mentioned three phases of a tumor population. 
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